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В статье приведены результаты числовых экспериментов по практической эффективности метода 
Брауна-Робинсон и монотонного алгоритма для матричных игр с точки зрения скорости сходимости и 
временных затрат, так и их точности. 
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У стати представлен! результати числових експериментв щодо практично! ефективност! метода Брауна- 
Робнсон 1 монотонного алгоритму для матричних 1тор з точки зору зы жности 1 часових обмежень та точност1. 
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Анализ конфликтных ситуаций приводит к возникновению новых моделей игро- 
вых задач. В работах [1-9] введен новый класс задач комбинаторной оптимизации 
игрового типа с двумя игроками, в которых платежи представлены матрицей, а на 
стратегии игроков накладываются комбинаторные ограничения. Такие задачи имеют 
широкое применение в разных областях человеческой деятельности. В [5] предло- 
жен приближенный итерационный метод решения таких задач, в основе которого 
лежат те же идеи, что и в основании метода Брауна-Робинсон [10]. Как известно, 
метод Брауна-Робинсон показал определенную практическую эффективность, но, на 
отдельных классах задач, как отмечено в [11] не достаточную скорость сходимости. 
Учитывая это для решения задач, рассмотренных в [1-9], возникает необходимость 
разработки других методов для решения матричных игр, которые будут основаны на 
других подходах. 

В [11] предложен монотонный метод решения игровых задач, который, по 
мнению авторов этой публикации, является более эффективным и имеет лучшую 
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скорость сходимости в сравнении с методом Брауна-Робинсон. В данной статье 
поставлена задача на основе числовых экспериментов сравнить известные методы 
решения матричных игр как с точки зрения скорости сходимости и временных 
затрат, так и их точности, чтобы использовать эту модификацию при решении 
задачи из [1-9]. 


Постановка задачи 


Рассмотрим классическую матричную игру [12], в которой первый игрок имеет 
т стратегий {= рт, второй игрок — и стратегий = рп. Игроки поочередно делают 
ходы: первый игрок выбирает свою 1-ю стратегию, второй — свою /-ю стратегию, 
после чего первый игрок получает выигрыш а, за счет второго игрока (если а, <0, 


то это означает, что первый игрок платит второму сумму |, ). 
Значение функции выигрыша можно записать в виде прямоугольной матрицы А: 


Я а Я 


а 


р они ЗЧ 
Таким образом, рассматривается матрица А=(а,), 1=1т, /=1Ъп размер- 
ности (т хп). Проведение каждой партии матричной игры с матрицей А сводится к 


выбору первым игроком 1-й строки, а вторым игроком /-го столбца, и получения 


первым игроком (за счет второго) выигрыша а,. 


Алгоритмы методов 


Если размерность исходной матрицы небольшая, то такую игру, как известно, 
можно свести к поиску седловой точки функции или к получению решения пары 
двойственных задач линейного программирования [12]. Матричные игры, в которых 
исходная матрица имеет довольно большие размерности, чаще всего решаются при- 
ближенными итерационными методами, в частности, методом Брауна-Робинсон [10], 
и монотонным итеративным алгоритмом [11] и др. Рассмотрим алгоритмы на- 
званных методов. 


Метод Брауна-Робинсон 


Идея метода — многоразовое фиктивное разыгрывание игры с заданной матри- 
цей А=(а,). Один розыгрыш игры называют итерацией, количество которых не- 


ограниченно. С ростом числа итераций М смешанные частоты применения чистых 
стратегией игроков приближаются к их смешанным оптимальным стратегиям. 

Расчеты делаются в предположении, что игроки хотят увеличить свой выигрыш 
(уменьшить проигрыш). Предполагается, что они не знают своих оптимальных 
стратегий. Ходы игроки делают в соответствии с принципом: будущее похоже на 
прошлое, учитывая все сделанные итерации. 

Первая стратегия (например, первым игроком) выбирается произвольно — возможно 
с целью увеличить возможный выигрыш. 
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Выбор следующего (5-го) хода 2-го игрока (после 5 ходов 1-го игрока) про- 


5 


5 к 
исходит выбором стратегии /, такой, что У‘, „=шшр) а, =у(5)-5 = 5, где 
` мы М4 та и 
(1, л),....(›./.) — пары стратегии игроков на /, ..., 5 ходах. А выбор (5 +1) -го хода 


1-го игрока, (после 5 ходов 2-го игрока) происходит выбором стратегии Г, согласно 
5 5 
условию: У`а_‚ =тахУ`а, =У(5):5=75. Цена игры у, как известно [10], удовлетво- 
5 /м 195т Ум 
М№М=1 М=1 
ряет неравенствуу(5) <у <у(5). 
Наглядно итерации метода изображены в табл. 1, где № — номер итерации, 
А, А,,...,А„ — соответственно стратегии первого, а В,,В.,...В, — второго игроков, 


Му ‚Му — минимальное в столбцах В, и максимальное в столбцах А, на шаге М; 


т 


в строку 2 на каждой итерации записывается накопление за все итерации (пла- 


Уи х 
тежи); у* = а — приближенное значение цены игры после № итераций. 


Окончание работы алгоритма может производиться при 1 выполнении равенст- 


ва Ух =Уты =У ; 2) после выполнения определенного количества итераций, в таком 


т. 


случае находят таху(М№) =у 


—_ тах 
1<М№<5 


и шшу(М№) =Ушь. Если У» <Ушь, ТО можно брать 
1<№<5 


У=(И к Увы) /2. 
Пример трех итераций решения игры с матрицей 
А 1 
А=|3 2 
14 
методом Брауна-Робинсон представлено в табл. 1 


Таблица 1 — Три итерации решения матричной игры методом Брауна-Робинсон 


мМ| 1 В | В || А | | 4 | Му у у* | му | у 
1 1 й) 1 о 1 и 4 1 1 2.5 4 4 
у 2 в [У Г 2 4 
2 3 1 4 1 2 3 1 3 1.5 2 5 2.5 
у 3 ЗЕ: 5 5 
2 3 2 1 2 3 1 6 2 2.3 8 2.7 
3 
у б о ПЕ: 8 б 


Рассмотрим другой итеративный процесс — монотонный итерационный алгоритм. 
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Монотонный итерационный алгоритм решения 
матричных игр 


Опишем согласно [11] монотонный алгоритм как итерационный процесс, кото- 
рый позволяет найти оптимальное значение у цены игры для игры /’,, которая за- 
дана матрицей А = (а, размерности тхпи. 

На нулевом шаге первый игрок выбирает произвольную чистую страте- 


2°т 


гию х' = (беж ря ) и определяется вспомогательный вектор с°как вектор скаляр- 
т ЕТС 
м 0 . 
ных произведений столбцов матрицы А и вектора х’, то есть с ‚ = У‘а,х, у =|п, 
1=1 


0 0 
с= (с С | 


п 
Шаг 1. Устанавливаем номер итерации М, равный 1, М =1. 


ы 5 ыы М „М1 . М-—1 
Шаг 2. Определяем у’ = шие” ти обозначим Л ={ а } множество 


У=Л, 


№-—1 м М 
индексов, на которых достигается у ‚тоесть “= Атешшс, . 
= 1<)<п 


Шаг 3. Формируем матричную игру Г” < Г, как частичную игру Г, с матри- 
цей А“ = (а), где 1=1т, ТЕ 
Шаг 4. Решаем сформированную частичную игру и определяем оптимальную 


=м =м м 
стратегию х =|^и,..., Хи |. 


=Мм =м =м =Мм и = 
Шаг 5. Определяем вектор с =(& абв следующим образом: с; м - 
151 
Ул=Ьи 
№-1 №-—1 

. Г ЗС. 

Шаг 6. Рассматриваем игру с матрицей размерности 2хп: | у у 
С1 ,...›Сп 


Находим оптимальную стратегию (а еб: ) ‚ О<а, <1 первого игрока в этой игре. 


При этом @ выбираем согласно следующему соотношению: 
и 3 м 
О: ат тпах ип [(1-ас ‘+ас; } 
а ый 


> М№М-—1 
Шаг 7. Если а, =0, то выход из алгоритма с ценой игры у , иначе — переход 
на следующий шаг алгоритма. 


М 
т 


=м 
М М М М1 
Шаг 8. Вычисляем значение вектора х’ =|[х,,....х |: хо =(1-@а,)х’ +аух . 
1 № М 


М № № 
Шаг 9. Определяем компоненты вектора с = вы ) согласно формуле: 


6 = (1- ау) с” +а,с”. 


Шаг 10. Увеличиваем номер М№ итерации на 1, переходим на шаг 2 алгоритма. 
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Числовые эксперименты 


С целью исследования практической эффективности рассмотренных методов для 
решения матричных игр и перспектив их применения в комбинаторной оптимизации 
для задач [1-9] они были программно реализованы на языке Иер. С использованием 
этой программной реализации методов проведена серия числовых экспериментов на 
компьютере с процессором АМРО ЕХ- 6300 с тактовой частотой 3.5 Ггц и объемом 
оперативной памяти 8 ГБ. 

Результаты экспериментов приведены в табл. 2 - 4, в которых указана раз- 
мерность задачи — параметры и и т - и время, которое было потрачено на 
вычисления (решалось по 10 задач при каждой размерности, приведено общее время 
вычислений всех задач), где «мс» — миллисекунды, «с» — секунды, «мин» — минуты, 
«ч» - часы. Элементы матрицы А для каждой задачи были сгенерированы как целые 
числа с равномерным распределением в интервале [1;1000]. Работа алгоритмов 
проводилась так: 

1. Для метода Брауна-Робинсон до достижения относительной точности по 


У ша 


цене игры вычислений = =10Ё—5 ‚т.е. при < = осуществлялась остановка. 


у тах 


2. Для монотонного алгоритма до получения значения параметра © =0, что 
позволило получать точное решение задач. Следует отметить, что монотонный алго- 
ритм требует решения на определенных этапах игровых задач размерности пх2, 
которые в данной программной реализации были решены с использованием симплекс- 
метода, что в некоторых случаях не позволило достичь максимальной эффектив- 
ности алгоритма. 

Для проверки точности решения были использованы решения этих же задач 
симплекс-методом, путем приведения игровых задач к паре двойственных задач ли- 
нейного программирования. 


Таблица 2 — Результаты числовых экспериментов (задачи типа их2) 


Размерность задач Время вычислений 
№ НЕ Монотонный 
пит п т АУ итерационный Симплекс-метод 
Робинсон 
алгоритм 
1. 1000 2 202 мс 68мс 16мс 
2. 5000 р 702 мс 302мс 22мс 
3. 10000 2 1с 631мсе 618мс 38мс 
4. 20000 2 2с 706мс 1с 231ме 63мс 
5. 30000 2 5 91мс 1с 897мс 86мс 
6. 40000 2 6с 657мс 2с 593мс 145мс 
7. 50000 2 6с 746мс Зе 237мс 178мс 
8. 100000 2 16с 569мс 6с 415мс 415мс 
9. 250000 2 42с 301мс 16с 672мс 1с 140мс 
10. 500000 2 1мин 3с 33се 851мс 2с 498мс 
11. 1000000 2 Змин 57с 1мин 30с 4с 924мс 
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Увеличение размерности задач в числовых экспериментах, до больших значений, 
чем те, что приведены в табл. 2 и 3, оказалось невозможным в связи с ограни- 
ченными вычислительными ресурсами компьютера (объемом оперативной памяти). 
В числовых экспериментах, результаты которых приведены в табл. 2, было решено 
ограничить размерность значениями и = т = 1000 в связи со значительным временем 
вычислений уже при этих параметрах. 

С помощью математического пакета СигуеЕхрен [13] была получена регрес- 
сионная зависимость времени работы алгоритмов от размерности задач и коэффициент 


2 
корреляции А“ (рис. 1-3). 
Зависимость времени вычислений от размерности(Мх2) 


250 
У= 2Е-10х2 + 6Е-05х + 3,3569,* 
В? =0,9916 „^^“ 


ю 

8 
` 

` 


® Браун-Робинсон в Монотонный « Симплекс-метод „ 


> 
| 
[=] 
` 
` 


100 де" У= 4Е-11%2 + 5Е-05х + 0,4234 , 
2 


вые ро оызибайй у = 6Е-1452 + 5Е-Обх - 0,0373 


> 

Е 

т 

= 

9 

® 

Ух 

= 

Е 

5 # 
5 

> 

о 

о 

[3 

3 

Е: 2 

В В? = 0,9994 
[= 


Размерность задач 
Рисунок 1 — Зависимость времени вычислений от размерности задачи 


(для случая их2) 


Таблица 3 — Результаты числовых экспериментов (задачи типа пхт) 


Размерность 5 
Время вычислений 
№ задач 
я п т Метод Брауна-Робинсон а Симплекс-метод 
итерационный алгоритм 

1. 10 10 1с 436мс 18 мс 12мс 

2. 20 20 7с 786мс 26мс 13мс 

3. 30 30 21с 446мс 95мс 24мс 

4. 40 40 38с 455мс 264мс 40мс 

5 50 50 1мин 29с 4с 317мс 84 мс 

6. 100 100 4мин 59с 2мин 47с 1с 123ме 
7. 200 200 20мин 5 с 16мин 20с 16с 182мс 
8. 300 300 42мин 24с 38мин 17с 1мин 25с 
9. 400 400 58мин 10с 1ч 5мин 4мин 10с 
10. 500 500 1ч 40мин 1ч 58 мин Эмин 50с 
11. 600 600 2ч 10мин Зч 40мин 24мин 10с 
12. 700 700 2ч 57мин 5ч 3Омин 42мин 20с 
13. 800 800 Зч 50 мин 7ч 3Омин 1ч 56мин 
14. 900 900 5ч 24мин 9ч 56мин 4ч 32мин 
15 1000 | 1000 6бч 59мин 18ч 20мин 7ч 8мин 
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Зависимость времени вычислений от размерности (МХМ) 


2Е 


70000 
60000 У= ЭЕ-05х? - 0,0547 + 19,313х - 701,68 " т 
В? = 0,9913 
5000 
Е ® Браун-Робинсон и Монотонный 4 Симплекс-метод 
© 40000 г 
> 
> = 
Е 
5 3000 ” 
> 
Е У ‚обл - 3,2538х +271,424 
Г: 5 В? = 0,9949 ле 2”” 
Е 20000 =” = =>“ 
| ВР И 
см „*' 
КО - У=9Е-05е: 5.079552 + 18,759х -597,77 
2 = 
о 


Размерность задач 


1000 


Рисунок 2 — Зависимость времени вычислений от размерности задачи 


(для случая ихт) 


Таблица 4 — Результаты числовых экспериментов (задачи типа 2хт) 


Размерность задач Время вычислений 
№ 5 
м я р Метод Брауна- Монотонный Симплекс- 
Робинсон итерационный алгоритм метод 
|. 2 100 892мс 35мс 19мс 
2 2 500 9с 634мс 582мс 244 мс 
3 2 750 13с 836мс 1с 408 мс 594мс 
4 2 1000 26с 492мс 11с 878мс 1с 309мс 
5 2 2000 1мин 4с 31с 109мс 5с 700мс 
6 2 3000 2мин 7с 57с 106мс 11с 721мс 
7 2 5000 4мин 33с 2мин 26с 36с 959мс 
8 2 10000 7мин 48с 5мин 44с Змин 2с 
9 2 15000 14мин 28с 14мин 8с бмин 34с 
10. 2 20000 З6мин 25с 29мин 11с 16мин 37с 
а Зависимость времени вычислений от размерности (2хМ) 
у= 5Е-10х3 - Ех + 0,0908х - 43, 959» 
ре = 0,9984 и 
у = 5Е-06х2 - 0,0128х +252 224 „ 
В? = 09977” Не" 


1500 


® Браун-Робинсон = Монотонный ^ Симплекс-метод, 


1000 


500 


Ей У 2Е-10% - 22:06 + 0,0153х - 11,332 
ый = 0,9974 


0 ОрЗЕН-х ЗИ МИАН 
0 


Время вычислений (секунды) 
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В результате проведенных числовых экспериментов выяснилось, что в случае раз- 
мерности задач 2хи и их2 решение монотонным алгоритмом проводится за меньшее 
время, чем методом Брауна-Робинсон. При этом следует заметить, что монотонный 
алгоритм, в отличие от метода Брауна-Робинсон, позволил получить для каждой за- 
дачи точное решение, но при этом определить вероятности стратегий только одного 
игрока (в методе Брауна-Робинсон — обоих). Для задач с разными значениями т и п 
при размерности, большей 400х400, вычисления монотонным алгоритмом прово- 
дятся дольше, чем методом Брауна-Робинсон. Это может объясняться необходимостью 
проводить часть вычислений с использованием симплекс-метода, что обусловлено 
данной программной реализацией. 

В ходе экспериментов выяснилось, что метод Брауна-Робинсон требует наименьший 
объем оперативной памяти, но при этом проводит значительно большее количество 
итераций. Также следует отметить, что в осуществленных числовых экспериментах 
решение задач симплекс-методом сведением матричной игры к паре двойственных 
задач линейного программирования, потребовало меньше всего времени, но в то же 
время использован значительный (наибольший из всех методов) объем оперативной 
памяти вычислительной системы. 


Выводы 


В статье приведены результаты числовых экспериментов по практической эф- 
фективности метода Брауна-Робинсон и монотонного алгоритма для матричных игр. 
В ходе проведенных исследований были проведены сравнения этих двух мегодов 
решения матричных игр С решением этих задач симплекс-методом с целью выясне- 
НИЯ ВОЗМОЖНОСТИ ИХ применения ДЛЯ решения задач комбинаторной оптимизации 
игрового типа. Полученные результаты позволяют сделать вывод, что монотонный 
алгоритм для задач большой размерности предпочтительнее, чем метод Брауна- 
Робинсон. 

Как направление дальнейших исследований можно рассматривать разработку 
методов решения задач комбинаторной оптимизации игрового типа с использова- 
нием идей, на которых основывается монотонный алгоритм. 
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О.О. 1етез5, О.М. ОШоузху, О.Т. ОШоу5Кауа 


Сотратбоп оГМешоа$ юг 5о[лиз Сате Та5к5: 


№ шитетса! Ехреттеп5 

ВасКогоипд: ш сотбтаюпа| орипиханоп а пе\ с1а5$ оР соттаюпа| орнпихайоп 
ргоетз оЁ Ше гатше фуре 15 сопз1Аегед, \уйй мо р/ауегз, ш уе раутепё аге 
ргезеще4 аз а таблх, ап4 сотЫтаюпа! гезбчсНопз аге ипрозе4 оп е зтаёеслез оЁ фе 
рауегз. Засв ргоепа$ сап Бе улаеу изе4 ш уаточ$ Не!4$ оЁР битлап асйуцу. То зоуе 
015 с1а5$ оЁ ргоептз ап арргохитае зо!уше тео4 1$ зиесехеа, \шсь 1$ Базе оп фе 
зате 14еаз аз Возе Бет? ш Фе Ба$15 оРе тефо4 оЁ Вго\т- Ко Мпзоп. И 15 \еП Кпо\п 
Фа Фе тефо4 ог Вгомп- Ко изоп Ваз зВо\уп а сейаш ргасйса1 еЁИслепсу, Во\уеуег ш 
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зерагайе с1аз5ез оЁ ргоетз сопуегоепсе зрее4 \’а$ пой зиЁЙслепте. Такшо 11$ шо сопз1- 
деганоп Юг зо1уше сотЬтаюпа! орипихайоп сате буре ий Бесотез песеззагу ю 4еуеор 
офег тео4$ Юг зо]уше тайлх сатез, уе \мШ Ъе Базе4 оп о ег арргоасВез. 


Маженча5 ап@ тео9д$: ш опг ага е уе Вауе зеё Ше 1азК ю сотраге Ше Кпо\п 
тефо4$ Юг зо[уше таблх гатез Базше оп патегса| ехрегипеп{$, Бой ш {егиаз оЁ зрее4 
ор сопуегеепсе ап4 {ите-сопзитриоп, аз уе аз ш {егил$ оРеиш ассигасу, ш ог4ег ю изе 
15 то сайоп а зо1уше сотбтаюна| орнпиханоп ргоепл$ оЁ хате фуре 


Вези! 65: ш ог4ег ю шуезисае Фе ргасйса! еЁесйуепез$ оР {Ве шейо@ оЁ Вго\т- 
Во п5$оп ап топоюшс а|еотйбта Юг 5оушФ таблх гатез аз ме аз Фе ргозресё5 Юг 
ег изе ш сотЫпаюпа[ орипихайоп ‚ Феу \еге паретете4 ш Ше 1апоиазе оЁ Оефы. 
Озше #15 зоЙ\аге ппрететайноп ог те#о4$5, а зетез оР питенса| ехрегитет Вауе 
Бееп ассотрИзВе4 оп а сотрщег УИ ап АМО ЕХ 6300 ргосеззог \/ИВ а с1осК зрее4 оЁ 
3.5 @НА апа 8 СВ оРКАМ. Ш Фе сопгзе оЁ Фе ехрегитеп$ 10 {а5К$ Вауе Бееп зо]уе4 ш 
еасб Апптеп$юп, Фе А таблх еетеп Юг еасВ {азК \уеге сепегайеЯ аз пцесег$ ий 
ип Ноги 9156бийоп ш Фе пиегуа[[1;1000]. То сВеск Фе ассигасу оЁ Фе зоИоп Фе 
знир!ех те#о4 Ваз Бееп изе4 фо зоуе Ше зате ргоеитз, Бу Бгиеше Фе гате ргоетз 
{0 а раш оЁ диа| ргоШепа$ ог Ппеаг ргозтапитте. УМН Фе Бер ог таетайса| расКасе 
СигуеЕхрее \е Вауе оМате4 геотеззлюп 4ереп4епсе оЁ Те уотК Ите оЁ а!еогИи$ оп 
фе Чипепзюп ое 1азКз ап4 1е сотге]аноп соеЁйслете А”, уме зВо\е4 а е ите 
Дерепа$ аз а зесоп4-4естее ро[упопта1. 


Сопсш$101: Ассогашо 1ю Ше гезиН$ оЁ ехрегипеп$ топоютс аеогИбт Юг Тагое- 
4пптеп$1оп ргоепа$ 1$ ргеегте4 1ю Фе тео4 оЁ Вго\уп- КоБзоп. 


Статья поступила в редакцию 26.12.2013. 
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